Spin state transition and covalent bonding in LaCoOa 
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We use the dynamical mean-field theory to study a p—d Hubbard Hamiltonian for LaCo03 derived 
from ab initio calculations in local density approximation (LDA+DMFT scheme). We address the 
origin of local moments observed above 100 K and discuss their attribution to a particular atomic 
multiplet in the presence of covalent Co-O bonding. We show that in solids such attribution, based 
on the single ion picture, is in general not possible. We explain when and how the single ion picture 
can be generalized to provide a useful approximation in solids. Our results demonstrate that the 
apparent magnitude of the local moment is not necessarily indicative of the underlying atomic 
multiplet. We conclude that the local moment behavior in LaCoOa arises from the high-spin state 
of Co and explain the precise meaning of this statement. 
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LaCoOa and related compounds have been much 
studied for half a century^ - — Their strongly tempera- 
ture dependent magnetic and transport properties have 
been eluding complete theoretical description so iai£-~— 
An apparent band insulator below 50 K, LaCoOa ex- 
hibits a local moment magnetic response above 100 K 
while the charge gap continuously disappears between 
450-600 KAir— This behavior points to an important 
role played by the electronic correlations as common 
among transition metal oxides. The correlated nature 
of LaCoOa reveals itself also in the formation of atomic- 
scale inhomogeneities with large magnetic moments upon 
moderate hole doping ! 14 ' 15 The picture of thermal evo- 
lution of LaCoOa as an entropy driven crossover from a 
non-magnetic to a magnetic state of Co 3+ ion has been 
generally accepted. However, the actual why's and how's 
are far from settled. The main open questions include: 
(i) Which atomic multiplet is responsible for the forma- 
tion of local moments? (ii) Does the lattice thermal ex- 
pansion actively contribute to the spin-state transition 
or is it merely a slave to the changes of the electronic 
structure? (iii) Why do the crossovers to local moment 
paramagnet and to bad metal happen at different tem- 
peratures? In this work we use the combination of the 
density functional band structures with the dynamical 
mean-field theory, known as LDA+DMFT, to address 
the former two questions. We discuss in detail the attri- 
bution of local moment behavior to a particular atomic 
multiplet in systems with covalent bonds, a question of 
general importance in oxide physics. 

The magnetic susceptibility of LaCo03 is usually an- 
alyzed in terms of the lowest multiplets of an isolated 
Co 3+ ion in octahedral crystal fields the low spin (LS) 
1 Ai (< 6 e°), the intermediate spin (IS) 3 Ti (iV) or the 
high spin (HS) 5 T 2 (t 4 e 2 ) states. The energy differences 
between the multiplets are controlled by the crystal-field 



(CF) splitting and the intra-atomic exchange J. While 
the LS singlet ground state is undisputed (at least for 
low temperature crystal structures) the nature of the first 
excited state is still a subject of debate. Goodenougbii 
attributed the appearance of local moment to popula- 
tion of the HS state. Heikes et alX on the other hand 
proposed IS scenario, which became popular— ~— after 
Korotin et alA obtained an IS ground state for expanded 
lattice with LDA+U calculation, contrary to a simple 
ligand-ficld theory. More recent experiments make a 
strong case for the HS scenario. The electron spin reso- 
nance shows a triplet excited state with a large g factor of 
3.35-3.55 which is consistently explained in the HS sce- 
nario invoking the effect of spin-orbit coupling The 
X-ray absorption spectra (XAS) and magnetic circular 
dichroism at the Z/2,3 edge of Co— also select the HS sce- 
nario. However, several author s 10 ' 22 " — pointed out that 
in order to interpret the magnetic susceptibility, specific 
heat, or XAS data in the HS scenario a rather strong 
temperature dependence of the crystal-field has to be as- 
sumed. In other words, the experimentally deduced in- 
crease of the HS population is significantly slower than 
anticipated with a fixed crystal field. The implication is 
that the apparent crystal field grows with temperature. 
This is somewhat unexpected since the expansion of the 
Co-0 bond which accompanies the spin crossover should 
reduce the crystal field. A possible explanation of this 
puzzling feature is provided by an inter-atomic repulsion 
between the HS states, which is equivalent to attraction 
between HS and LS states . 10 ' 23 ' 24 Breathing lattice dis- 
tortion proposed by Raccah and Goodcnough, 3 studied 
in detail by Bari and Sivardiere^ 5 - provides a mechanism. 
Knrzek et al— used LDA+U calculations to argue for 
HS-LS attraction. Recently a purely electronic mecha- 
nism of HS-LS attraction was observed by two of us— in 
a two band Hubbard model. 
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It is well known that the effective crystal field in tran- 
sition metal oxides is largely due to hybridization with 
ligands, i.e. generated by hopping of predominantly e g 
electrons between the metal and oxygen sites. This effect 
is particularly pronounced in LaCoC>3 and leads to charge 
fluctuations on the Co site. Therefore it seems natural 
to question descriptions based on an isolated ion. How 
to define and distinguish HS and IS states when sizable 
charge (valence) fluctuations are present? Or more pre- 
cisely, is it possible to express the T-dependent suscep- 
tibility as a sum of contributions from different atomic 
states? We will demonstrate that in a solid the answer 
is negative in general. However, in the insulating phase 
of LaCo03 the notion of LS, IS and HS states can be 
preserved when these are generalized to include the hy- 
bridization induced fluctuations. The T-dependent para- 
magnetic moment can be to a good accuracy approxi- 
mated by a T -independent HS moment multiplied by a 
T-dependent weight. Importantly, the magnitude of this 
effective moment differs from the free ion value. 



The role of the lattice thermal expansion poses a 
chicken-and-egg question about the relationship between 
the spin-state transition and the anomalous lattice ex- 
pansion i 20 ' 28 ! 29 While the transition to the magnetic 
state, both IS and HS, weakens the Co-0 bond and thus 
leads to its expansion, stretching the Co-0 bond reduces 
the effective CF splitting and thus favors a magnetic 
state. Therefore there is a positive feedback between 
these two effects. To include the lattice response directly 
into our calculations is not computationally feasible at 
the moment. Therefore we have performed calculations 
for several lattices corresponding to experimental crystal 
structures at different temperatures. By comparing the 
T-dcpendence of the spin susceptibilities on these lat- 
tices we find that experimentally observed variation of 
Co-0 bond-lengths has pronounced effect on the elec- 
tronic properties. 



Inclusion of thermal effects in 'first principles' density 
functional approaches is notoriously difficult. Therefore 
the existing studies are either limited to the T = LS 
phased or assume that temperature enters only through 
the lattice thermal expansion4i2£ The only serious at- 
tempt to explicitly include temperature in such a calcu- 
lation was made by Edei— using the variational cluster 
approximation (VCA). The VCA and DMFT methods 
share many formal similarities. However, while DMFT 
allows calculation of the one-particle as well as multi- 
particle correlation functions, VCA is limited to the one- 
particle correlations. As pointed out in Ref. [Io|, the VCA 
spin susceptibilities or atomic state occupancies are ob- 
tained from a reference C0O6 cluster. The relevance of 
these quantities for an infinite system is a conjecture, 
which calls for verification with other methods. 



I. METHODS 
A. Model 

Multi-band Hubbard Hamiltonian with the two- 
particle interaction within the Co 3<i shell is used to 
describe LaCo03. The one-particle part of the Hamil- 
tonian, which spans the Co 3c? and O 2p orbitals, has 
been constructed from the local density approximation 
(LDA) to the density functional theory. The non-spin- 
polarized bandstructure obtained with WIEN2k2i was 
transformed into the Wannier basis representation with 
wicn2wannier— and wannier90^ codes. The Hamilto- 
nian in this representation reads 

H = ( h k d a pdl aa d k ^ + h^ s pl 1<7 p kS a+ 

h(T 

(i) 

Here dk aa and Pfc 7<T are Fourier transforms of the an- 
nihilation operators di aa , pi^ a which destroy the d or p 
electron with the orbital index a or 7 and the spin index 
a in the ith unit cell, nf aa . are the corresponding occupa- 
tion number operators, and h k :b aj3 are the corresponding 
matrix elements of the one-particle LDA Hamiltonian. 

The U°p are the density-density matrix elements in 
the e g -t 2g basis of the full Coulomb interaction param- 
eterized with U (Slater parameter Fq) and Hund's ex- 
change J. The screened values of U = 6.0 eV and 
J = 0.8 eV have been obtained using the constrained 
density functional theory (DFT) scheme.— ~— Calcula- 
tions, where U and J were varied, are also reported to 
assess the stability of the results. 

The h dd diagonal elements were modified to account for 
the static part of the interaction, double-counting correc- 
tion, 

4% = ftjgj - (N orb - 1 ) UnS aS , (2) 

where n is the average self-consistent occupancy per Co- 
rf orbital, U is the orbital averaged interaction energy, 
and N or i, is the total number of interacting orbitals on 
single site (10 in our case).— This is equivalent to sub- 
tracting the orbitally averaged Hartree potential felt by 
the d electrons. 



B. DMFT selfconsistency and one- particle spectra 

The one-particle Green's function of the Hamiltonian 
(H|) is found by iteratively solving the DMFT equations on 
the Matsubara contour. The auxiliary impurity problem 
is solved by the continuous time quantum Monte-Carlo 
(CT-QMC) method in the hybridization expansion for- 
mulation^ using the implementation based on free-access 
package ALPS i 39 ' 40 The Wang-Landau re- weightin g 41 ' 42 
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was employed in order to ensure the ergodicity of the 
simulations for some parameter values, in particular at 
low temperatures and close to the spin state transition. 

Once the calculation was converged we have evaluated 
the one-particle spectra in real frequency and analyzed 
the impurity dynamics and spin susceptibility. For the 
former analytic continuation is necessary To this end 
we have employed the maximum entropy method in two 
modes, (i) continuation of the local Green's function from 
the imaginary time r to real frequency uj, (ii) contin- 
uation of the local self-energy from the Matsubara fre- 
quency iuj n to lu. For the latter we have used the statis- 
tical error estimates following Ref. H^. The former was 
used to crosscheck the results of (ii) and the spectra are 
not shown in the paper. Analytic continuation of self- 
energy has several attractive features, such as being ex- 
act in the non- interacting limit, providing a direct access 
to the fc-resolved and ligand spectra, and smearing out 
the lifetimes (imaginary part of the self-energy) but not 
the quasiparticle dispersions. 



C. Susceptibility and local moments 

To analyze the on-site dynamics we have studied two 
additional quantities. First, the state weights^ which 
are the diagonal terms of the site-reduced density matrix, 
i.e. expectation values of the projection operators on the 
atomic (many-body) states = \/j,)(fi\. In the present 
case of the density-density interaction the site-reduced 
density matrix is diagonal in the occupation number ba- 
sis. Therefore the knowledge of the state weights allows 
us evaluate expectation value of any local operator, e.g. 
the instantaneous local moment (ml) = m 2 z (p)(P il ) 
with m z ([i) = (n\rh z \(i). 

Second, we define the imaginary time state-state cor- 
relation matrix C^^r) and its time average LT^,, 



CW(r) = ft(r)P„(0)) 

n M „ =T / drC^(r), 
Jo 



(3) 



where (3 = 1/T is the inverse temperature. The cor- 
relation matrix allows us to analyze the local response 
functions (via fluctuation-dissipation theorem) and de- 
compose them into different contributions. In particular, 
we will be interested in the local spin susceptibility \, 
which in the paramagnetic state can be expressed as 



X= dT(m z (T)m z (Q)) 
Jo 



(4) 



The second expression shows that in general the local 
magnetic response cannot be decomposed into contribu- 
tions of atomic states, but pairs of states must be consid- 
ered. Decomposition into states contributions is possible 



only in a special case, such as an isolated atom where 
C m „(t) can be made diagonal. 

Discussing briefly the physical meaning of these quan- 
tities we start by pointing out that in the course of time a 
given atom visits various quantum-mechanical states as a 
result of statistical (thermal) fluctuations and quantum- 
mechanical (causal) evolution. The weight of a given 
state is a relative measure of the time spent by the atom 
in this state, which does not distinguish between thermal 
fluctuations and causal evolution. The state-state cor- 
relations distinguish to some extent between these two 
effects as only states connected by causal evolution can 
have a non-zero cross term. The state weights can be 
obtained as row (or column) sums over iI M „. 



II. RESULTS AND DISCUSSION 

A. Non-interacting band structure 

We have considered the experimental distorted per- 
ovskite structure with R3c space group containing two 
formula units per unit cell. The structural parameters 
were taken from the X-ray measurements of Ref. [l8|. To 
assess the effect of lattice thermal expansion the calcu- 
lations were repeated for the experimental structural pa- 
rameters obtained at three different temperatures (de- 
noted as Tiattice in the following): 5, 450, and 750 K. 
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FIG. 1: The calculated orbitally resolved spectral function 
[0:p (shaded cyan), Co:tg 9 and e 3 (red and black line)]. Left 
panel: the hybridization term h!^^ of the Hamiltonian (TJJ) 
was set to zero. Right panel: the LDA Hamiltonian in the 
Wannier basis. 

The octahedral crystal field splits Co:<i states into six 
t2 g states at lower energy and four e g states at higher 
energy. The resulting orbitally resolved spectral density 
is shown in Fig. [TJ The splitting is strongly contributed 
by the 0:p-Go:d hybridization, as is clear from the com- 
parison in Fig. [TJ where the right and left panels show 
the spectral function with and without the p-d hybridiza- 
tion included. The p-d hybridization increases the dis- 
tance between the centers of ti g and e g bands consider- 
ably. The band broadening, more pronounced for the e g 
band, is another consequence. Matching 0:p and Co:<i 
features reflect formation of bonding and anti-bonding 
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states. Stronger hybridization of the e g orbitals com- 
pared to t2 g ones results from their larger spatial overlap 
with O.p orbitals. The ti g band is further split to the 
doublet and the a\ g singlet due to a distortion from 
the octahedral symmetry. This splitting does not play 
an important role in our study, though. As in previous 
calculations the metallic ground state is incorrectly pre- 
dicted by LDA£ 



B. Thermal effects and lattice expansion 

LDA+DMFT calculations were performed for several 
temperatures between 290 and 2320 K (/3 from 40 to 
5 cV _1 ). If not stated otherwise, the results are shown 
for U = 6.0 eV and J = 0.8 cV (obtained by cLDA cal- 
culations). 

Local susceptibility. The local spin susceptibility Xi 
calculated from Eq. 21 is shown in Fig. as a func- 
tion of T. For T[ attice = 750 K the huge error bar at 
T = 580 K is due to a long autocorrelation time despite 
the Wang-Landa u 41 ' 42 sampling. For all lattice parame- 
ters we observe an emergence of Curie-like susceptibility 
at high temperatures. The corresponding average spin 
moments m scr = \JT\ are shown in Fig. [5Jd. As we cal- 
culate directly the local susceptibility (the local response 
to a field applied to a single site of the infinite crystal) 
the inter-site exchange does not enter the definition of 
the local moment. Fig. [2b suggests a gradual thermal 
population of a magnetic state. The lattice expansion 
clearly favors the magnetic state and the experimentally 
observed magnitude of the Co-0 bond-length expansion 
has a sizable impact on our results. 

Next, we analyze x(T) for a fixed lattice. Assuming 
the isolated atom scenario the local susceptibility can be 
written as 
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where [i and v are the magnitude of the Ising moment 
and the multiplicity of the excited magnetic state, A e fi 
is its energy with respect to the LS ground state. In 
Fig. [2J: we show A e g derived from the calculated suscep- 
tibility using v = 6 and (i = 3.5, which correspond to an 
Ising HS state, a choice explained later in the text. Like 
in the experiments 2 ^— and the VCA theory^ our A e g 
for a fixed lattice increases with temperature for reason 
discussed below. 

It is quite clear that our results do not provide a good 
quantitative description of LaCoOa as the spin crossover 
takes place at a too high temperature. This is not surpris- 
ing since the present theory is unlikely to achieve the nec- 
essary ~10 meV accuracy without fine tuning the mate- 
rial parameters by hand. However, two important trends 
are revealed. First, the lattice response (expansion of 06 
octahedra around moment-carrying sites) acts as a pos- 
itive feedback for generation of local moments. Second, 
this is countered by a purely electronic effect, making 



FIG. 2: (a) The local spin susceptibility as a function of 
the temperature for different lattices ( Tiauice)- (b) The cor- 
responding screened spin moment. The dashed line indicates 
the value for the IS in the ionic limit. The letters A, B, C 
denotes the solutions discussed in the text, (c) The apparent 
crystal-field A e g obtained from Eq. [5] 



addition of a local moment the harder the more local 
moments are already in the system; this is reflected in 
the increase of A c ff with the temperature (Fig. [2}:) . This 
is another way of saying that there is a repulsive inter- 
action between the magnetic sites in the LS background. 
We point out that in our calculation all Oq octahedra are 
the same, which excludes a possible contribution of the 
breathing distortion. 3 Indeed, Kyomcn et al.— substitut- 
ing Co with Al, Ga and Rh came to the conclusion that 
electronegativity rather than ionic radius of the neighbors 
is the parameter which controls A e g. We suggest the fol- 
lowing picture based on the observation that a strong 
Co:e s -0 bond favors LS state, and that in the Co-O-Co 
trimer the Co-0 bonds share the central p a orbital. Due 
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to this sharing the energy gain per Co-0 bond in the 
trimer is less than the energy gain for an isolated Co-0 
bond. Therefore breaking (weakening) one bond in the 
trimer makes the other bond stronger. Introducing a lo- 
cal moment on one Co site provides this bond breaking 
and strengthening the other bond favors the LS state. 



eating holes to dominate the electronic transport. Alike 
the VCA results^ the photoemission part of the spectra 
exhibits a transfer of spectral weight from the low-energy 
peak (~1 eV) to higher energies (~3 eV) observed exper- 
imentally^ 
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FIG. 3: (Upper panels) The orbital-resolved spectral func- 
tions [states/eV/formula unit], t% g (e 9 ) : solid red line, tg 9 
(dig): dashed red line, e g : black line, 0:p: cyan shaded area. 
The 0:p spectral function is downscaled by the factor of 2 to 
fit in the graph. (Lower panels) fc-resolved spectral function 
A^(lo) along the high symmetry directions. Left panels show 
the interacting non-magnetic (low-spin) solution (denoted as 
A in Fig. [2b; T = 580 K, Ti atUce = 5 K) and the right pan- 
els display the paramagnetic solution (with a large content of 
high-spin atomic states) (denoted as C in Fig. [2b; T — 1160 K, 

TiatUce = 750 K) . 

Spectral functions. In Fig. [3] we compare the one- 
particle spectra of the low-T non-magnetic state (left 
panels) and the high-T paramagnetic state (right). The 
orbital resolved spectra are displayed in the upper panels 
and the fc-resolved spectra along the high symmetry di- 
rections are shown in the lower panels. The non-magnetic 
spectrum resembles the LDA solution (Fig. [TJ, the main 
difference being a uniform (Hartrec) shift of the e g band. 
There is very little dynamical renormalization since the 
LS state is an approximate cigenstatc of both the kinetic 
term (dominated by the crystal field) and the interaction 
term taken separately. The correlated nature of LaCoC>3 
is revealed at elevated temperature. The thermal popula- 
tion of the excited atomic multiplets leads to a formation 
of local moments, which are incompatible with disper- 
sive bands. As a result incoherent features appear in the 
spectra. The nature of the charge gap changes from a 
semiconductor like gap between coherent valence ti g and 
conduction e g bands (left panel of Fig. [3]) to a ti g -ti g gap 
(right panel of Fig. [3]). The bottom of the conduction 
band is now defined by an incoherent t2 g excitation, the 
tail of which gradually fills the gap with the increasing 
temperature. The top of the valence manifold is formed 
by a renormalized dispersive ti g band. This spectrum is 
consistent with the positive Seebeck coefficien t 45 ! 46 indi- 



C. Spin state analysis 

Local state statistics. Next, we address the local dy- 
namics at Co sites and whether it is meaningful to de- 
scribe it in terms of the atomic states (such as LS, IS, 
HS). The hybridization expansion CT-QMC solver is well 
suited to this task as it provides the site-reduced statisti- 
cal operator (density matrix), referred to as state statis- 
tics^ This quantity describes the probability of finding 
an atom in a particular many-body state and the expec- 
tation value of any local operator can be easily obtained 
from it. We display schematically the atomic states im- 
portant for the forthcoming discussion in Fig. 01 
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FIG. 4: Atomic states at LS and HS configurations. The 
cyan lines depict the p-oxygen orbitals and the black cir- 
cles denote the hole in p-O shell. The effects of the crystal 
field and Co:d-0:p hybridization are schematically depicted. 
In total four atomic states belong to LS block (see the left 
part of Figure): (LS) all ttg orbitals occupied, all e g orbitals 
empty, (LS+e 9 ) in addition to LS single e g orbital occupied, 
[LS+2e 9 (t4-)] in addition to LS two e g orbitals with the op- 
posite spin orientation occupied, [LS+2e 9 (tt)l i n addition to 
LS two e g orbitals with the parallel spin orientation occu- 
pied. All these states are accessible from LS via the e g hy- 
bridization. For brevity we do not distinguish LS+2e 9 (ti) 
from LS+2e 9 (tt) from now on. In total three atomic states 
belong to HS block (see the right part of Figure): (HS) all 
majority spin orbitals and one minority spin t% g occupied, 
(HS+e 9 ) in addition to HS another minority spin e 9 orbital 
occupied, (HS+tg 9 ) in addition to HS another minority spin 
ttg orbital occupied. All these states are accessible from HS 



There are many atomic states with non-negligible 
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weights contributing to the partition function (Fig. [5|. 
The total contributions of different charge states (insets 
of Fig. [5]) point to sizable valence fluctuations, which is 
related to finding substantial admixtures of d 7 L and d 8 L 2 
state to the d 6 ground state in the cluster calculations^ 
Unlike the cluster calculations, in DMFT the ligand hole 
does not remain coherent with the central Co atom due 
to the influence of the rest of the crystal. Therefore one 
cannot use the CoC>6 eigenstates to analyze the local dy- 
namics. Instead, we use the statistical description and 
also analyze the temporal evolution of the atomic states. 

Clearly, the d 6 atomic multiplcts denoted as LS, IS, 
and HS in Fig. [5] are not sufficient to describe the local 
physics in LaCo03. We distribute the atomic states into 
LS, IS, and HS blocks (Fig. 2]). Although an a priori as- 
signment of the blocks is not unique we later present an 
a posteriori justification of our choice. For example, the 
state denoted as HS+^g can, in principle, be reached by 
adding a t2 g electron to the HS state as well as by adding 
an e g electron to the IS state. In Fig. [5] we present the 
state statistics for various lattice parameters and tem- 
peratures corresponding to the solutions denoted as A, 
B and C in Fig. [5b. Besides substantial weights of the 
d 7 states the figure reveals that the increasing local mo- 
ment response is connected to the growing weight of the 
HS block, while the IS block has only minor weight. 

Correlation matrix of local states. Although we have 
identified the atomic states with large weights, a question 
arises whether their appearance is due to a unitary evolu- 
tion or rather due to statistical averaging. This question 
is closely connected to the decomposition of susceptibility 
into atomic states contributions fEqs. [3l4)l . States con- 
nected by a unitary evolution lead to a large off-diagonal 
element of the time averaged state-state correlation ma- 
trix n and thus their individual contributions to the sus- 
ceptibility cannot be well defined. On the other hand, if 
the weights of different states originate in statistical aver- 
aging the corresponding off-diagonal element of n is van- 
ishing as is its contribution to the susceptibility. Analysis 
of the correlation matrix n^,,, displayed in Figs. [5] and [7] 
reveals non-zero off-diagonal elements, indicating a uni- 
tary evolution between the corresponding states. Nev- 
ertheless, the matrices can be arranged in a block diag- 
onal form, which justifies our choice of the LS, IS and 
HS blocks. As the unitary evolution between the most 
populated LS and HS blocks has a vanishing probability, 
their simultaneous population is a result of the statisti- 
cal averaging and blocks generalize the notion of atomic 
multiplets in isolated atoms. 

The block-summed contributions to the local suscep- 
tibility are indicated by numbers inside the respective 
blocks in Figs. [51 [7] In both cases the contribution of 
the HS diagonal block amounts around 97 % of the total 
susceptibility. Inspecting the block contributions to the 
spin-spin correlation (Fig. [8]) we find a finite r-constant 
part of the HS contribution, leading to Curie-type sus- 
ceptibility, in contrast to the rapidly decaying LS contri- 
bution. This allows us to define an effective HS moment 
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FIG. 5: The weight of dominant atomic states for U = 
6.0 eV, J = 0.8 eV. (a) non-magnetic solution (A in Fig. [2b, 
T = 580 K, r lattlce = 5K), (b) low-T solution (B in Fig. [2b, 
T = 580 K, T latttce = 450 K), (c) high-T solution (C in Fig. [2b, 
T = 1160 K, T lamce = 750 K). 
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We obtain /zhs of 3.52 and 3.56 fiB in the low-T and the 
hight-T solutions respectively. The weak T dependence 
of the effective moment and its dominant contribution to 
the susceptibility x justifies expressing x as a product 
of Curie term /i 2 /T and a T-dcpendent weight. Cova- 
lent Co-0 bonding results in about 10% reduction of the 
effective moment from its atomic value of 4 [Lb- 
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FIG. 6: The correlation matrix 11^ between the dominant 
atomic states for T — 580 K and Ti attice = 450 K (solution 
B of Fig. [2]d). Color-coded values show the state- by-state 
relative contributions in % to the sum over all pairs. The 
numbers within blocks of atomic states indicate the contri- 
bution of the blocks to the local susceptibility in the units of 
10" 3 emu • Oe" 1 ■ mol" 1 . 
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FIG. 7: The same as Fig. © for T — 1160 K and n 
750 K (solution C of Fig. [3>). 
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D. Role of U, J, double counting correction 

Since the form and construction of the Hamiltonian 
(1) is to some extent an ad hoc procedure it is impor- 
tant to understand the sensitivity of our conclusions to 
the particular values of U, J, and the double counting en- 
ergy (Edc is the second term on r.h.s. of Eq.[2]). Although 
these are not unadjustable parameters in the present the- 
ory, their determination is not unique, which holds in par- 
ticular for the double counting correction E dc . In Fig. [§] 
we show the dependence of the weight of the LS block 
for various values of U and J. As expected J is the more 
important parameter. Its cLDA value falls into the spin 
state crossover range of 0.8-0.9 eV. Outside this range 
the results are insensitive to temperature or the variation 
of U. Variation of U inside the crossover regime has some 
impact as higher U suppresses the fluctuations to the d 7 L 
and d 8 L 2 states. As a result the LS state is destabilized. 

We have examined the role of double counting correc- 




l - 
0, 



x(eV ) 



FIG. 8: The HS (upper panel) and LS (lower panel) nor 
malized block contributions to the spin-spin correlation func 
tion, J2 
1160 K and 

^lattice 



750 K (solution C of Fig.fjp). 




FIG. 9: Weight of the LS block of atomic states (LS, LS + e 9 , 
LS + 2e 9 ) for various U and J at Tiattice =5K, (a)T= 580 K, 
(b) T = 1160 K. 



tion (E dc ) for U = 6.0 eV, J = 0.8 eV, T = 1160 K, 
nattice = 750 K. We varied E dc by SE dc in the range 
of ±5eV around the self-consistent value of 33.56 eV. 
The positive values of 5E dc mean that the Co:d states 
are shifted downwards in energy closer to the 0:p state, 
which in turn enhances the hybridization and leads to 
preference of a LS metallic phase. The weights of the 
LS, IS and HS blocks are shown in Fig. [TU] For SE dc > 1 
the system becomes metallic and the definition of the LS, 
IS, and HS blocks loses its justification. 

Based on the one-particle spectra and the overall be- 
havior of our results we conclude that the self-consistent 
value of E dc provides a rather good description of the ac- 
tual material. We also point out the T-dcpendcnt vari- 
ation of E dc is rather small and had minor effect on the 
T-dcpcndcncc of both \ and A e ff. 



III. CONCLUSION 

We have studied the temperature dependence of the 
magnetic and spectral properties of LaCoOs using the 
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FIG. 10: The weight of LS, IS, and HS blocks of atomic 
multiplets for various double counting energies for high-T so- 
lution (U = 6eV, J = 0.8 eV, T = 1160 K, r !ttttice = 750 K). 
The values indicate the offset from the self-consistent value of 
33.56 eV. 



LDA+DMFT approach. Our results show that the local 
moment response at elevated temperatures is associated 
with the HS state of Co ion and that there is an effective 
inter-atomic repulsion between HS atoms in the LS ma- 
trix. Our findings agree with the VCA calculations^ and 
LDA+U cluster calculations^ 6 - as well as with the conclu- 
sions of the experimental studies^ - — On the other hand, 
our results are inconsistent with interpretation of the in- 
termediate temperature phase in terms of IS state d 18 ! 20 
Furthermore, since a purely electronic mechanism of the 
HS-LS attraction exists^ the absence or smallness of 
the breathing mode distortion^ does not exclude HS- 
LS short range ordering. The experimentally observed 
(anomalous) lattice expansion has a pronounced effect 
on the LS-HS crossover, which leads to the conclusion 
that the lattice provides a strong positive feedback for 
the crossover. 

To account for the strong covalent bonding with O the 
notion of Co HS state has to be generalized to include 
not only d 6 , but also d 7 and d 8 electron configurations. 



This leads to reduction of the magnitude of effective HS 
moment. Therefore using the apparent magnitude of 
the magnetic moment as an indicator of the underlying 
atomic state may lead to incorrect conclusions. 

Although we have varied the computational parame- 
ters in a wide range we have not found a phase that 
would be dominated by the IS state. Therefore scenar- 
ios invoking HS-IS crossove r 24 ! 26 are not consistent with 
our results. Guided by a similar LDA+DMFT study on 
metallic SrCo03, it is plausible that in the metallic phase 
of LaCo03 observed experimentally above 600 K distinc- 
tion between IS and HS state is not possible as those are 
connected by unitary evolution of the system. 

Based on the above observations we propose follow- 
ing scenario of LaCo Oa ph ysics, which is in many as- 
pects similar to Refs. Il0ll26l . At the lowest temperatures 
LaCo03 is LS with isolated Co ions in the excited HS 
state. Increasing temperature assisted by the lattice feed- 
back leads to growing density of the HS sites. Effective 
repulsion keeps the HS sites apart leading to a short range 
HS-LS order, which is responsible for the insulating be- 
havior in the 100-500 K range, similar to the observation 
made in Ref. [27]. The second crossover around 500 K is 
associated with 'melting' of the LS-HS order. This leads 
to an anomalous lattice expansion due to the breaking 
of attractive LS-HS bonds. The experimentally observed 
onset of mctallicity changes the local moment character 
by coherently admixing IS-like states to the dominant HS 
configuration. The distinction between HS and IS in the 
hight T metallic phase is thus not possible. 
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